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ABSTRACT 

In supernova spectroscopy relatively little attention has been given to the properties 
of optically thick spectral lines in epochs following the photosphere's recession. Most 
treatments and analyses of post-photospheric optical spectra of supernovae assume that 
forbidden-line emission comprises most if not all spectral features. However, evidence 
exists which suggests that some spectra exhibit line profiles formed via optically thick 
resonance-scattering even months or years after the supernova explosion. To explore 
this possibility we present a geometrical approach to supernova spectrum formation 
based on the "Elementary Supernova" model, wherein we investigate the characteris- 
tics of resonance-scattering in optically thick lines while replacing the photosphere with 
a transparent central core emitting non-blackbody continuum radiation, akin to the 
optical continuum provided by decaying 56 Co formed during the explosion. We develop 
the mathematical framework necessary for solving the radiative transfer equation under 
these conditions, and calculate spectra for both isolated and blended lines. Our compar- 
isons with analogous results from the Elementary Supernova code SYNOW reveal several 
marked differences in line formation. Most notably, resonance lines in these conditions 
form P Cygni-like profiles, but the emission peaks and absorption troughs shift redward 
and blueward, respectively, from the line's rest wavelength by a significant amount, 
despite the spherically symmetric distribution of the line optical depth in the ejecta. 
These properties and others that we find in this work could lead to misidentification of 
lines or misattribution of properties of line-forming material at post-photospheric times 
in supernova optical spectra. 
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Introduction 



Many of the physical processes which contribute to spectrum formation in a supernova (SN) 
change dramatically as it ages. Given the difficulty and complexity of including all such processes, 
as well a s their evolution in time , exploration into this computational frontier has begun only very 



recently (I Pinto k Eastman 



2009; De et al 



2010; 



2000 



Hoflich 



Dessart k Hillier 



201 



2003 

T 



Kasen ct al 



2011 



2006 



Jack et al. 2009; Kromer fe Sim 



Hillier k Dessartll2012l ). Most work in SN spec- 



troscopy has focused on the early "photospheric" phase, comprising the time from explosion to 



SDectrum (Branc 


l et al. 


19831. 1 


58,4 .Tefferv et al. 19921: Mazzali k Lucv 


1993: 


Mazzali et al. 


1995|; 


Fisher et al. 


1997 


Mazzali et al. 


1997: 


Hatano et al.l 


1999al: 


MazzalilboOO. 


2001 


Branch et al. 


2005: 


Dessart k Hillier 


2005; 


Mazzali et al. 


2005; Hachineer et al. 2006). Consider ablv less attention 



has been given to the "nebular" phase, several months or more after explosion, where emis- 



sion f rom forbidden lines constit ut e most of the spect r um ([ Axelrod 



1992 



Mazzali et al 



2011 



Kozma k Fransson 1998al lbl; iMazzali et al. 1998 ; Maeda et al 



1980; Ruiz -Lapuente fe Lucyl 



2006; iJerkstrand et al 



2011 



Finally, the inter mediate phases, th at is, a few months post-explosion, have 

One reason so much attention focuses 



received the least amount of scrutiny (jMaurer et al. 



2011 



on either very early or very late times is that in these regimes one can reproduce with reasonable 
fidelity the physical processes which dominate spectrum formation through a number of si mplify- 
ing a ssumptions. For example, in the early, photos pheric phase the Soboley approximation (jCastor 
1970j) and a resonance-scattering source function (jjefferv fe Branch! Il990l ) are both accurate ap- 
proximations due to the high densities and steep velocity gradients in the SN ejecta; at late times 
one typically assumes bot h that the ejec ta is optically thin and that line emission arises exclu- 
sively from forbidden lines (|Axelrodlll980l ). However these two groups of assumptions are generally 
incompatible with each other in the intermediate regime of a SN. 

Complicating matters further is the possibility that the evolution of the different types of 
line-forming processes is asynchronous. Specifically, the emergence of forbidden emission lines in a 
SN spectrum may not herald the systematic withdrawal of resonance-scatteri ng in optically thick 



permitted lines. Though this possibility has frequently been acknowledged (jBowers et al. 



Branch et al. 



2005 



Jerkstrand et al 



2011 



1997 



it has to our knowledge never been pursued in detail 
until now. To illustrate this point we show in Figures Q] and [2] SYN0W fits to the day +87 optical 
spectrum of SN 2003du and the day +115 spectrum of SN 1994D, respectively. In the SN 2003du 
fit we use only three ions: Na I, Ca II and Fe II. The photospheric velocity is 7000 km s , and the 
photospheric temperature is 8000 K. The excitation temperatures of all three ions is 7000 K. In 
the SN 1994D fit, we include five ions: Na I, Ca II, Cr II, Fe II, and Co II. Here the photospheric 
velocity is 6000 km s -1 and the temperature is 10,000 K. The excitation temperature of all ions is 
7000 K. The observe d and synthet i c spe ctra of both SN 1994D and SN 2003du have been flattened 
using the method of IJefferv et al.l (|2007l ). 



The fits to the observed spectra in both of these figures are relatively good, and in both fits 
only permitted lines were considered. Even among the features which SYN0W cannot reproduce 
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accurately, most of the emission peaks and absorption troughs in the synthetic spectra form at the 
same wavelengths as in the observed spectra, and only the strengths of the features are disparate. 
(The exceptions to this are the wavelength range 6600 - 7800 A in SN 2003du and 6600 - 7500 A in 
SN 1994D, throughout which SYN0W fails completely to reproduce the observed features.) It would 
therefore be a remar kable coincidence if the observed features in these two objects arise from purely 
forbidden emission ((Branch et al.ll2005l ). 



In addition to the circumstantial evidence for persistent resonance-scattering provided by these 
reasonably good SYN OW fi ts to several-months-ol d type la supernovae (SN e la), calculations by 
Hatano et al.l (|1999bl ) and Ijerkstrand et al.l (|201ll ) further defend this claim. lHatano et al.l (| 1999b! ) 
calculated the Sobolev optical depth (Equation 1 of that paper) for the most common ions observed 
in SN la optical spectra. In Figure 9c of that work they show the Sobolev LTE optical depth for Fe 
I & II and Co I & II for an iron-rich composition at 80 days post-explosion in a SN la model, and 
find that the optical depths for those ions can be as high as r ~ 50, 000. Similarl y, the simplified 
form of the Sobolev optical depth shown in Equation 39 of Ijerkstrand et al.l (120111 ) shows that the 
optical depths for some UV resonance lines (e.g., Mg II A2798) can still be as high as r ~ 1000 
even eight years post-explosion. 

In light of the evidence presented above for resonance-scattering at relatively late times in 
SNe la, we present a ge ometrical explorat i on of this phenomenon in the spirit of the Elementary 
Supernova approach of iJeffery &; Branch! yj)90). Rather than attempt to discern exactly how 
late after a SN explosion resonance-scattering continues to contribute significantly to spectrum 
formation, we assume simply that the photosphere has receded and that the continuum opacity 
in the core of the SN is negligibly small. We replace the photosphere with a transparent core 
emitting non-blackbody continuum radiation and distribute ions both inside this "glowing" core 
and outside the core in a non-emitting shell. We then calculate emergent flux profiles for lines with 
parameterized optical depths in several different geometrical configurations. We scrutinize line 
formation and behavior in both isolated and blended cases. Because SYNOW inspired this model, 
the two naturally invite comparison, which we indulge throughout this work. In $2] we treat the 
case of pure continuum (no lines) in the core; in $3] we introduce a single line in the core; in §@] 
we surround the glowing core with a transparent, non-glowing shell and move the line to the shell; 
in <|5] we account for the effects of multiple scattering by placing two lines in the core; $6] contains 
our discussion of the implications of our results on interpretation of late-time spectra of SNe; in 
£J7] we conclude our work. Finally, we include Appendices which contain complete mathematical 
derivations so that all of our work may easily be reproduced. 



2. A transparent, uniformly emitting core 

We begin by assuming the SN is spherically symmetric. Under this circumstance it is natural 
to work in (p, z) coordinates, where p is the impact parameter of a ray relative to the center of 
the SN, such that the p = ray exactly bisects the SN; and z is orthogonal to p, with the z = 
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Fig. 1. — A SYNOW fit to the +87 day spectrum of SN 2003du. The synthetic spectrum contains 
Na I, Ca II, and Fe II, each with an excitation temperature of 7000 K. The photospheric velocity 
is 7000 km s 1 , and the spectrum has been divided through by the blackbody continuum. The 
permitted line identifications (in blue) are from the SYNOW fit, with Fe II fea tures unlabeled. Th e 
alternative forbidden line IDs (in green) follow those in other SNe la made by I Bowers et al.l (|1997l ). 




Fig. 2. — A SYNOW fit to the +115 day spectrum of SN 1994D. The synthetic spectrum contains 
Na I, Ca II, Cr II, Fe II, and Co II, each with an excitation temperature of 7000 K. The photospheric 
velocity is 6000 km s , and the spectrum has been divided through by the blackbody continuum. 
As in Figure Q] the unlabeled features in the SYNOW spectrum are due to Fe II. 
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line also exactly bisecting the SN. We will at times transform to spherical polar coordinates for 
computational expedience, where the radial coordinate r satisfies 

2 2,2 ri\ 

r = p + z . (1) 

We also assume that the SN undergoes homologous expansion, v = r/t, so that surfaces of constant 
line-of-sight velocity are vertical planes, that is, planes of constant z. We further assume that the 
observer is located at z — > — oo, in which case all rays incident on the observer are parallel. Next, 
we assume that the photosphere has receded and has been replaced by a transparent, spherically 
symmetric core with outer radius r c which emits continuum (but not blackbody) radiation. The 
post-photosphe ric spectra of many S Ne la contain a flat and weak continuum which is either thermal 
in nature or, as 



Bowers et al.l (|1997l ) suggest, due to the sea of weak optical lines of lowly ionized Fe 
and Co. At very late times when the spectrum becomes truly nebular, there is no optical thermal 
continuum, but we do not seek to extend our methods into this very late regime. To mimic this 
pseudo-continuum we assign to the core a grey, spatially constant volume emissivity j\(r) = j for 
all radii r <r c and all wavelengths A. This assumption allows us to present simple analytic results; 
with only slight modification, our methods for calculating line profiles are amenable to chromatic 
emissivity. 



2.1. Continuum only 

In general, to calculate the emergent flux from a SN atmosphere one must first calculate the 
source function at all locations in the SN, followed by the emergent intensity of rays exiting the 
ejecta toward the observer. However, since in this model we neglect all continuum opacity, the 
source function is not a well defined quantity in the absence of lines. Therefore we write down 
immediately the emergent intensity of a constant p ray originating at the bacl>0 of the core and 
traversing toward the observer, without calculating the source function. When no lines are present 
the intensity of a ray passing through the core is proportional to its geometric length: 

I{p)=2j{rl-p 2 ) l l\ (2) 

A representative intensity ray is shown in FigureO If instead the emissivity has spatial dependence, 
j = j( r )> the result has a more general form: 

lip) = / j((p 2 + z 2 ) l ' 2 )dz = 2 / j((p 2 + z 2 f/ 2 )dz (3) 



( r 2_ p 2)l/2 J Q 

where we have transformed the argument of j since the intensity along a ray of constant p depends 
only on its z coordinate. 



In this discussion "front" and "back" refer to locations in the SN nearest to and farthest from the observer, 
respectively. 
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Already a physical peculiarity arises: limb darkening in the absence of scattering or absorption. 
The traditional interpretation of limb darkening in photospheric objects such as dwarf stars and 
young SNe is that intensity rays with large impact parameter p, that is, rays which emerge from 
the limb, reach an optical depth of r ~ 1 at shallower layers in the atmosphere than low-impact 
parameter rays. Under the assumption of local thermodynamic equilibrium (LTE), shallower loca- 
tions in atmospheres have lower temperatures and thus weaker source functions, since LTE requires 
by definition that S = B where B is the Planck function. A weaker source function in turn leads 
to a lower intensity emergent ray, the cumulative result of which is limb darkening. 

In the case of this post-photospheric model, however, high-p rays accumulate less continuum 
as they proceed from the back of the object toward the observer. Such accumulation does not occur 
in photospheric-phase atmospheres of SNe except through lines in the atmosphere which scatter 
continuum photons from the blackbody photosphere into the ray; however, this contribution is 
small except when the line has an extraordinarily high optical depth. 

Using Equation [2j the continuum flux at all wavelengths is given by 

f 2TT f r ' c 47T 

F x = h cos 6dn = ~y / hpdp = —jr c . (4) 

J r c JO 6 

For clarity, and without loss of generality, hereafter we set j = 1, so I\, J\, and Fx have units of 
length and the geometric nature of our results becomes apparent. In addition, since in homologous 
expansion v oc r, we interchange lengths and velocities where convenient. 



2.2. Mean intensity 

We now calculate the mean intensity, J, both inside and outside the glowing core. Without 
any lines, J is not a particularly useful quantity in this model since we do not need it to calculate 
the emergent flux. However, when we add a single line to the core in <|3] we require the source 
function to calculate the emergent intensity, and in the resonance-scattering approximation S = J. 

By definition, 

Mr) = ^J h(r)dQ = \j UrW (5) 

where [i = cosO, and we have applied to the radiation field the condition of azimuthal symmetry. 
Inside the core, r <r c , so from Figure 

X = rfi + (rV + r 2 c - r 2 ) 1/2 . (6) 

Setting I\(r) = X(r), plugging this into Equation [5j and using standard integral tables, we find 

J( r ) = — \rr c + ^ -± In 



r c + r 



(7) 
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Calculating the mean intensity outside the core (with the line still inside the core) is slightly more 
complicated and we include the derivation in Appendix [Aj The result is: 



rr c + 



(r 2 



In 



r + r c 



(8) 



Figure [5] s hows the behavior of J(r) in units of r c and for comparison the dilution factor W(r) 
()Mihalas!ll978l ) is also shown. The shape of J is vary similar to that of a Gaussian and it is larger 
than W(r) until quite large r, where both functions behave as \{^) 2 • We emphasize that J peaks 
strongly at r = 0, a result which affects line formation significantly, as we will discuss in $31 



3. A line in the core 



We now treat the case of a single line, forming inside the core. We assume, for simplicity, a 
constant Sobolev optical depth r in the line throughout the core. Before continuing, we note that 
the most profound effect of homologously expanding matter on the radiation field is that photons 
redshift with respect to the matter regardless of propagation direction. Therefore photons in a ray 
which starts at the back of the core and emerges toward the observer redshift continuously as they 
move along the ray. Referring to Figure (U if a ray originates in front of the plane of constant z res 
corresponding to the location in the core z = z res where the line Doppler shifts into resonance with 



a particular wavelength point — that is, if p > 



a/2 



then the intensity of that ray is 



simply its continuum value, given by Equation [2j However, if the ray forms behind that plane 



if p < 



a/2 



the line attenuates some of the continuum intensity by scattering photons 



out of the ray and therefore out of the observer's line of sight. In the Sobolev approximation this 
attenuation manifests as a e~ T term multiplying the continuum intensity at the location of the line. 
In addition to attenuating the intensity along a ray, the line also contributes to the intensity via 
its source function S(r[z res \). Specifically, the contribution is 



I line = 5(1 - e" T ). 

The two different types of rays are depicted in Figure El Their emergent intensities are 



(9) 



2(rl 



.p2)l/2. 
. ^1/2 



3 )e- T + S(r)(l-e- T ) + (r 2 



P 



2x1/2 



+ Z r . 



p < (rl 
p > irl 



a/2 



2 

res/ 

2 a/2 

res/ 



(10) 



The integral over I Zrcs {p) is straightforward except for the term containing the source function S{r 
To integrate this term we define 

f "\/ r c~ Z Tcs f'*\/ r c ^res 

2r 



g[Zres = 



S(r)pdp 



rr c + 



In 



r c + r 



pdp. 



(11) 



We then transform the integration variable from p to r, which leads to 



1 

2^ 



rr c + 



In 



r c + r 



rdr. 
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Fig. 4. — Geometric configuration used to calculate J(r < r c ) both for pure continuum and for a 
single line in the core. X is the magnitude of the vector ab, r c is that of a&, and r is that of be. 




Fig. 6.— 



Intensity rays with a single line in the core. 
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This integral is still unwieldy, so we transform integration variables once again by denning 

■2-res 



fi = 



and 



from which we find 



Mr 



' Mres 

r c M 

We then change the integration variable from r to fi, which yields 



fibres) = -wl4 



1 + 



M 



Mres 



2[j,[i r 



In 



The result is 



„3 r 



g(Zves) = -f 

b 



1 + 2 In 2 



Mres 



Mres /o 



Mr 



In 



1 + Mr 
1 - Mr 



1 + Mr 
1 - Mr 



dfj, 



ln(l 



Mres/ 



(12) 



The remainder of Equation H] poses little challenge and leads directly to an analytic result for the 
flux profile for the constant-r case: 



2tt 

2 



— + (1 

3 v 



") 



g(Zres) + 



s) 



|3l 1 



(13) 



We remark here that, although Equation [J] defines the flux as a function of wavelength, we have 
derived Equation [13] in terms of the coordinate z res . One transforms between z o A using the 
first-order Doppler formula, 

A-A 



An 



(14) 



where r max (f max ) is the maximum radius (velocity) of the eiecta. We show in Figure [7] the line 
profiles of Na I D, A5892, with several different optical depthso For comparison we show in FigureE] 
the profile of Na I D as calculated by SYNOW, with the same optical depths as in Figure [71 

The profiles in both plots appear quite P Cygni-like, but differences do exist between them. 
We highlight two in particular. First, when a photosphere is present as in SYNOW, the emission peak 
of a line profile is located at its rest wavelength, regardless of the velocity of the material forming 
the line in the SN atmosphere. This fact s implifies enormously the task of identifying isolated lines 
in a SN spectrum (jJeffery &: Branchlll990l ). In the post-photospheric model, however, the emission 
peak is not centered on the line rest wavelength; rather it is considerably redshifted. As Figure [7J 
shows, with a core outer edge velocity of only 3000 km s~ • L , the emission peak is ~ 15 A redward of 



2 Examining Figures [T] and [2l identification of Na I is ubiquitous; we therefore use the Na I D line for illustrative 
purposes. 
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Fig. 7. — Flux profiles for Na I D, A5892, in the post-photospheric model, with spatially constant 
optical depth and an outer core velocity of 3000 km s -1 . The vertical dashed line indicates the 
rest wavelength of the line. The r = 10 and r = 100 profiles overlap almost exactly and are 
indistinguishable in this figure. In contrast to SYNOW the red emission peak does not occur at the 
rest wavelength, but rather is redshifted by an amount independent of optical depth. 



the line rest wavelength. Though the absorption component of a P Cygni profile is used more often 
than the corresponding emiss ion to identify lines in a spectrum (since the former dominates the 



latter in overlapping lines; see lJefferv &: Branch! 1 19901 ) . this redshift could lead to misidentification 



of lines in intermediate- or late-time SN spectra. 

What is the origin of the emission peak redshift? This question is easier to answer in the limit 
r — > oo, the line profile for which would look identical to that of r = 100 in Figure [7J In the 
r — > oo case the line at the resonance point z Tes scatters all continuum photons forming in the ray 
segment z > z Tes out of the ray, and the contribution of the line to the intensity along that ray 
is simply Iij ne = S(r[z Tes \). Thus only two sources of emission contribute to the flux at a given 
wavelength: 1.) the continuum emission in front of the resonance plane at z = z ies ; and 2.) the 
source function at the plane. If the plane is located near the back of the core, z < r c , the portion of 
the core's volume which emits continuum photons which escape unscattered and reach the observer 
is large. However, because the surface area of the plane is small when z is close to r c , and because 
from Figure we see that S is small at r ~ r c — only half of its maximum value — the emission 
contribution from the source function on the plane is in turn relatively small. 

As the resonance plane moves forward (to bluer wavelengths, smaller z), the volume of the 
emitting core in front of the plane decreases, but the surface area of the plane grows, and Figure [5] 
shows that the source function at the plane grows quickly as well. From ^r c < z < r c the scattering 
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Fig. 8. — Flux profiles for Na I D in SYNOW, with photo-spheric velocity 3000 km s and spatially 
constant optical depth. The maximum velocity of the line-forming region is 6000 km s _1 . The right 
dashed vertical line indicates the rest wavelength of the Na I D line (A5892), and the left vertical 
line indicates the blueshifted rest wavelength of Na I D at the photospheric velocity. The r = 10 
and r = 100 profiles overlap almost exactly and are indistinguishable in this figure. Unlike the 
profiles in Figure the emission peak in the SYNOW case remains fixed at the line rest wavelength 
regardless of the optical depth. 

emission from the resonance plane more than compensates for the diminishing continuum emission 
from the core, causing the flux to increase monotonically as one moves blueward in that region, 
eventually reaching the emission peak, which in Figure[7]is ~ 5915 A. Blueward of this peak, despite 
the fact that S increases monotonically until reaching z = 0, the surface area of the plane increases 
only slightly when z > and the now very large resonance plane occults so much of the core over 
a small shift Az that it can no longer compensate for the large amount of emission removed from 
the volume of the core behind it, causing the flux to decrease as z becomes smaller. This transition 
point exists redward of z = and thus the emission peak of the P Cygni profile is redder than the 
rest wavelength of the line. 

Blueward of z = 0, that is, — r c < z < 0, the area of the resonance plane begins to decrease, now 
only obscuring a cylindrical volume of the glowing core. Even though the length of this occulted 
cylinder increases as z becomes more negative, its radius decreases, and the portion of the total 
volume of the core that this cylinder comprises decreases as well. The now-unocculted limbs of the 
core, emitting continuum which the resonance plane can no longer scatter away, grow in volume 
and eventually compensate once again for the cylindrical volume obscured by the plane, causing 
the flux to increase. This transition point manifests as the absorption minimum in the flux profile 
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(~ 5855 A in Figure [7]), blueward of which the flux increases until the resonance plane reaches 
z = —r c and we recover pure continuum. 

The second difference between the line profiles with and without a photosphere is that, in the 
former, the flux in strong lines deviates from the continuum by a large amount, whereas in the 
latter the changes are small. The flux in the absorption minimum of the To = 100 line in Figure El 
for example, is 80% lower than the continuum value. In our post-photospheric model, on the other 
hand, even the strongest line in Figure [7] departs by only up to 30% from the continuum. That 
the post-photospheric model exhibits such small departures from continuum is due to the ability 
of a small portion of the core to emit unscattered photons toward the observer even when the 
resonance plane is near the front. In the SYNOW case, the shell above the photosphere does not emit 
any continuum, so when the plane is close to the front of the ejecta, only scattering from the plane 
itself contributes to the flux, causing it to be extremely low at the absorption minimum. 

The physical simplicity of constant r line profiles such as those in Figures [7J and [8] readily 
facilitates analysis of line formation in SN ejecta, as we have just seen. However, one must also 
be aware of complicating effects such as the inhomogeneous structure of SNe, viz. the velocity- 
dependent density profile. In the Elementary Supernova framework one assumes implicitly that 
r = r(p), where p is the matter density, and accounts for this by writing r as, e.g., an exponential 
or power law function which decreases with increasing velocity. To illustrate how a variable optical 
depth affects line formation we show in Figure [9] the line profile of Na I D in our post-photospheric 
model with t(v) = tq exp(— v/v e ), where tq is a constant and v e = 500 km s _1 . Again, for 
comparison, we show in Figure [10] a SYNOW plot of the same line with the same exponential r 
profile. 

In the post-photospheric case, shown in Figure [9] the relationship between velocities and the 
positions of the emission maximum and absorption minimum is not immediately obvious. Like the 
SYNOW line profiles, shown in Figure [TO] the blue edge of the absorption component (not of the 
absorption minimum; see 5835 A in Figure and 5775 A in Figure [TO]) indicates the maximum ve- 
locity of the line-forming material, although in blended cases this can be difficult to discern. Unlike 
SYNOW, in this model both the peak and trough continuously redshift and blueshift, respectively, 
as the line optical depth increases. In the SYNOW case, on the other hand, the separation between 
the rest wavelength of the line and the minimum of the blueshifted absorption component indicates 
the velocity of the photosphere and thus the minimum velocity of the line-forming materi al itself 



(assuming the material is not "detached" above the photosphere. iJefferv fe Branch! Il990l ). Only 
when the line optical depth is extremely high does the location of the absorption minimum begin 
to move blueward of the photospheric blueshift, as in the To = 100 line in Figure [TO] For each 
of the three weaker lines in that figure, one can simply measure the blue edge of the absorption 
minimum, ~ 5830 A, then calculate the photospheric velocity, ~ 3000 km s _1 . 

We also note that the r = 10 and r = 100 line profiles in Figures [7J and [8] appear identical, 
whereas in Figures [9] and [10] the tq = 10 and tq = 100 profiles are distinct. This difference is due to 
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Fig. 9. — Flux profiles for Na I D, A5892, in the post-photospheric model with outer edge core 
velocity 3000 km s _1 and a line optical depth following an exponential decay, t(v) = tq ex.p(—v/v e ), 
where v e = 500 km s _1 . The vertical dashed line indicates the rest wavelength of the line. In 
contrast to SYN0W, the red emission peak does not occur at the rest wavelength, but is rather 
significantly redshifted by an amount which depends on the strength of the line. Also unlike SYN0W, 
the blue absorption minimum blueshifts continuously with increasing optical depth. 

the velocity dependence of r: in the former pair of figures, the optical depth is constant everywhere, 
leading to a saturated profile throughout the line-forming region when r > 10. In the latter pair, 
on the other hand, r decreases exponentially with increasing velocity, so the line samples such large 
optical depths at only very low velocities. This velocity dependence affects line formation outside 
the core in exactly the same way (cf. $U). 



4. A line outside the core 

Having explored in detail the geometric effects of continuum emission on line formation in the 
core, we now introduce a non-emitting, transparent shell around the core, with outer radius r s , 
where r s > r c . In this region there is no continuum emission and line formation occurs in exactly 
the same way as in SYN0W. Therefore, unlike the case discussed in the line forming region now 
exists outside the core, rather than inside. We include the shell in this model to account for the 
possibility of intermediate mass elements such as Ca II forming lines at late times above the material 
which has been burned all the way to the iron-peak. 
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Fig. 10. — Flux profiles for Na I D, A5892, in SYNDW, with photospheric velocity 3000 km s -1 and 
a line optical depth following an exponential decay, t(v) = Toex.p(—v/v e ), where v e = 500 km s _1 . 
The maximum velocity of the line- forming region is 6000 km s _1 . The right vertical dashed line 
indicates the rest wavelength of the line, and the left dashed line indicates the wavelength of Na I 
D blueshifted to the photospheric velocity. Except for extremely high optical depths, the blue edge 
of the absorption minimum remains fixed at the photospheric velocity; and for all optical depths, 
the red emission peak remains at the rest wavelength of the line. 

4.1. Intensity for a line outside the core 

We must now consider a number of possible ways that the line and a particular ray can interact. 
First, if the ray has impact parameter p > r c then it never intersects the core and, regardless of 
the location of the line resonance point z res , the ray's emergent intensity is 

I Zl Jp) = S(r[z vcs ])(l-e' T ) (15) 

where S(r) is now the source function outside the core, given by Equation [8j and 

r = r(r[z res ]). (16) 

If the line has p < r c but z Tes < — \/r 2 — p 2 es , then the line attenuates the continuum ray from the 
core: 

IztJp) = S(rfeJ)(l - O + 2(r c 2 - p 2 ) 1 /^. (17) 
Finally, if the line has p < r c and z res > \/r 2 — p 2 es , then it does not attenuate the core continuum 
and the emergent intensity is 

h xe M = S(rW)(l " z~ T ) + 2(r c 2 -P 2 ) 1/2 - (18) 
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These three cases are shown in Figure [TT1 



4.2. Flux for a line outside the core 

While F\ proceeds simply from I\ in the case r < 77, its form is much more complicated in 
the r > r c regime. In particular the flux integral takes a unique form in five different regions. We 
show the five zones in Figure [12j 



1. —r 8 < z < —(rj - r 2 ) 
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Figure [131 depicts single- line profiles with different values of r. We first note the resemblance 
of these line profiles to those of SYNOW shown in Figure This is not surprising, given that the only 
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Fig. 12.— 



The five flux regions in the core+shell configuration. 
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difference between the two models is that in SYNOW the core is opaque while in the post-photospheric 
model is it transparent. Both exhibit strong deviations from the continuum flux value, especially 
in the absorption component. We discussed in $3]the cause of this near-zero flux in the absorption 
component of the line. 

There is one major difference between the line profiles in Figures [8] and [T3j the flat-topped 
emission component in the latter. To discern the source of this plateau in the spectrum, we conduct 
a geometric analysis similar to that presented in again studying the limit r — > oo. Starting with 
the resonance plane at the back of the shell, z = r s , we see that as it moves forward toward the 
observer, its surface area grows, which leads to the flux increase from 6010 A blueward to 5950 A 
in Figure [131 

When the resonance plane reaches the back edge of the core and begins to move forward 
through it, that is, when < z Tes < r c , the core projects onto the plane a central circular region 
where r = and continuum emission from the core is unattenuated. The component of the plane 
which samples the optical depth in the shell is a ring with area 



From this equation we see that A is constant when < z res < r c ; this constancy is the cause of the 
emission plateau in Figure [13l 

When — r c < z res < 0, the plane begins to obscure the core, starting at the core's limb, and 
scatters an increasing amount of the continuum emission out of the observer's line of sight. When 
Zres < — r c the entire core is obscured and the only emission from the ejecta comes from the source 
function at z res . In the region — r s < z res < — r c , line formation occurs in exactly the same way 
as in SYNOW, leading to the absorption trough in Figure [13] which is almost identical to that of the 
SYNOW profile shown in Fig. [8] 

The constant-flux emission in the specta in Figure [13] is the most distinguishing feature of 
our post-photospheric model. In ^6] we will consider a nebular line-forming region with a "hole" 
devoid of line optical depth, and there we will encounter a geometric conspiracy similar to the one 
presented in this section, leading to similarly flat features in line profiles. We remark in addition 
that SYNOW can produce flat-topped emission features in a spectrum by detaching lines from the 
photosphere. Thus we conclude that in general a plateau-shaped line emission feature indicates 
some kind of missing line opacity. 



We now study the case of two resonance lines in the core. Understanding the ways in which 
multiple lines "interact" via their overlapping flux profiles, and in turn being able to identify them 
individually, is critical to interpreting observed spectra since, in all epochs of a SN, its optical 
spectrum exhibits severe line blending. 
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5. 



Two lines in the core 
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Fig. 13. — Flux for Na I D in the shell, with core velocity 3000 km s and outer edge shell velocity 
6000 km s _1 . The rest wavelength of the line is Ao = 5892 A, indicated by the left vertical dashed 
line. The right dashed line indicates the rest wavelength of Na I D redshifted to the back edge of 
the core, at which point the emission peak forms a plateau. For r c < r < r s , r = constant and for 
t < r c , t = 0. The r = 10 and r = 100 profiles overlap almost exactly and are indistinguishable in 
this figure. 
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In the resonance-scattering approximation, for any number of lines in the core, the source 
function of the bluest line always has the single-line form given in Equation [7J since there are no 
bluer photons in the radiation field which can redshift into resonance with it. The second bluest 
line then interacts only with its bluer neighbor; the third bluest line interacts with its two bluer 
neighbors; etc. In computational terms this means that one calculates the source function for each 
line starting with the bluest and moving redward. The details of the calculation are shown in 
Appendix [Bj 

We first consider the case v COTe = 3000 km s" 1 , and we hold the optical depths of both lines 
fixed everywhere in the core, tb = tr = 1. We then set the rest wavelength of the blue line to 
2850 A and decrement the rest wavelength of the red line from 3000 A, where the lines are too far 
apart in wavelength space to blend in the spectrum, to 2855 A, where the lines overlap completely. 
The result is shown in Figure [TH In this figure we find that l ine blending in the post -photospheric 



case occurs in essentially the same way as in the SYN0W case (jJefferv &i Branch! 1 199d ). Specifically, 
the absorption component of the red line completely overwhelms the emission of the blue line, and 
also, when the two lines overlap perfectly (the bottom-right panel of Figure [EJ), the redshifted 
emission component appears the same as the single-line case, but the absorption is much deeper. 

One peculiar feature of Figure [14] is that, when the lines overlap only partially, the emission 
component of the blue line acquires a triangular shape. In fact, in the upper-right panel of that 
figure, the triangular emission peak is nearly symmetrical. However the abrupt (perhaps unrealistic) 
drop in r at the edge of the core in this calculation may somewhat exaggerate this triangular feature. 

In this section we have restricted our discussion to the interaction between two lines confined 
to the core. However, with the introduction of multiple lines, many other types of interactions are 
possible. As an example, in Appendix IB.2I we construct the geometric framework for the two- line 
case where both lines are outside the core. However the situation is more complex due to the myriad 
combinations of line locations as well as their combined interaction with the core and shell. We 
thus defer a thorough examination of all the possible effects that can occur in this case, including 
the effects of the photosphere before it has receded completely, to a future work (B. Friesen & 
R. C. Thomas, in preparation). 



6. Discussion 



Much recent work analyzing SNe la spe c tra has focused on t h e importance of asymmetries 



(Maeda et al. 



2010 



2011 



Maund et al. 



2010 



Foley Kasen 



2011 



Foley et al.1 12012 ). However 



the spectra produced by this model illustrate a potentially significant complication in interpreting 
spectra at late times in SNe la, specifically in the interpr etation of line emission in the absence of 



scattering. To illustrate this point we refer to the work of lMaeda et al.l (|2010l . l201ll ). who interpret 



observed blueshifts and redshifts of Fe lines in nebular spectra as resulting from asymmetric clumps 
of iron-peak material moving toward and away from the observer, respectively, in the optically thin 
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Fig. 14. — Interaction of two lines in the core as they blend together due to Doppler broadening. 
The rest wavelength of the blue line is fixed at Xb = 2850 A, while the rest wavelength of the red 
line moves blueward. The optical depths of both lines are fixed at tb = tr= 1. 
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SN la ejecta. This interpretation implies that the explosion of a white dwarf resulting in a SN la 
is itself highly asymmetric, but relies heavily upon the assumption that spectral features at late 
times arise only through forbidden transitions of atomic lines due to the low gas density, and that 
the resulting emission profiles are distributed roughly symmetrically about the rest wavelength. 
Indeed, this assumption has become widespread in nebular spectrum synthesis since the pioneering 
work of I Axelrodl ()1980l ). However, if an appreciable amount of resonance-scattering is present at 
this epoch, our results show that in the presence of significant continuum emission, a perfectly 
spherically symmetric distribution of matter produces a redshifted emission compone nt. It remains 



unclea r what effect, if any, that resonance-scattering has on the sample considered in lMaeda et al 



{[2010j) since the objects in their study were often several hundred days older than the objects we 
attempted to fit using SYNDW in Figures [T] and [2j 



Maund et al.l 1)20101 ) studied polarization measurements of the o bjects inlMaeda et al.1 ()2010l ) 



for which such data were available, and drew similar conclusions as iMaeda et al 



2010) , that is, 



that a symmetry alone may explain the diversity observed in SNe la. The survey in iMaund et al 
(|2010l ) was biased in that polarization data existed only for SNe with redshifted emission peaks for 
particular lines of Fe II and Ni II, but nevertheless we note that, at intermediate times, for many 
features there are multiple possible line identifications, and redshifted emission can result from 
either optically thin, receding ejecta, or, as shown in Figure [TJ, from spherically symmetric e jecta 
with large line optical depth. While the very late data (jMaeda et alj|2010l : IMaund et al.ll2010i ) may 
indeed be showing the inferred asymmetry, our results show that one should indeed be cautious 
about the epoch of the data when interpreting redshifted emission profiles. 

Finally, we consider another significant effect of resonance-scattering at late times in SNe. Stud 



iesof late-time line profiles and the imaging of the supernova remnant S Andromeda (jHoflich et al 



2004 



Motohara et al 



2006 



Fesen et al 



2007 



Gerardv et alj|2007l : IMaeda et alj|20ld ) have lead to 



the suggestion that a "nickel hole" exists in at least some SNe la. In optically thin media, a central 
region devoid of line-forming materia l manifests in a spectrum as a flat-topped emission feature 
in the spectrum; iHoflich et al.1 (J2004J) find just this when studying the [Fe II] A16, 440 line in a 
spectrum of SN 2003du taken ~ 300 days after explosion, suggesting that 56 Fe and therefore 56 Ni 
were absent in the central part of the SN ejecta. However, the IR observations in these studies are 
quite noisy, and if the same effect could be observed earlier, or in stronger lines, there would be 
more flux, allowing statistics to be built up on the existence of a nickel hole in SNe la, which would 
provide important constraints on the underlying explosion model. 

With this in mind we explore the possibility that one may observe a similar plateau feature 
inside the glowing core by carving out regions of zero optical depth in a resonance-scattering line. 
To study this, we return to the core-only model, with velocity 3000 km s _1 , and set the rest 
wavelength of a single line in the core to 3000 A. We then set the optical depth of the line to r = 1, 
and proceed to exclude this optical depth from progressively larger portions of the core's central 
region. The resulting line profiles are shown in Figure [Pol The classical rounded profile is replaced 
by a flat profile, but there still exists both an emission peak and an absorption trough. This result 
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corroborates the plateau effects we studied in £J5] — Figure [13] already showed that when a region 
of large line optical depth surrounds a region of zero optical depth, flat-topped components can 
appear in the spectral line profiles. 



7. Conclusions 

We have presented the geometric framework for a post-photospheric model of a SN, making 
simple assumptions about its emissivity and source function, in an attempt to explore the effects 
of resonance-scattering in optically thick lines in SN atmospheres in epochs where such scattering 
processes are rarely considered. Our model is inspired by and similar to the Elementary Supernova 
model, but the substitution of a photosphere for a transparent but continuum-emitting core leads to 
significant differences in line formation. The most noticeable difference is that the emission peak of 
a line is redshifted from its rest wavelength, even though the model is spherically symmetric. This 
property may affect the interpretation of SN spectra in post-photospheric phases, when asymmetric 
effects on line formation become influential. 

Our model remains mostly within the geometric confines of the Elementary Supernova model, 
and choices about its properties different than the ones made here are possible. For example, one 
may assume that severe line blending in the core creates a pseudo-continuum, leading to a source 
function in the shell exactly equal to that assumed in SYNOW: 

S = W{r)J, 

where W(r) is the dilution factor. We have also chosen to parameterize all line optical depths rather 
than calculating them in any self-consistent manner, e.g., by solving rate equations. We believe 
our choices in these matters have resulted in a level of detail commensurate with the simplicity of 
this model. 

To explore fully the other implications of resonance-scattering on line formation at post- 
photospheric tim es in SNe la we are currently incorporating the formalism presented in this work 



into SYNAPPS (IThomas et al.ll201ll ). whereby we will be able to analyze the interaction of the 
receding photosphere with the emerging glowing core of 56 Co and its effects on the SN spectrum, 
a process which has so far received little attention. 
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Fig. 15. — A post-photospheric line profile where the line is excluded from the various parts of the 
inner region of the core. The core has outer velocity 3000 km s _1 . The optical depth in the region 
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A. S(r) outside the core 



In 2] we calculate sample spectra for lines forming outside the glowing core. Therefore, we 
calculate the source function in the region r > r c . From Figure [TBI 



ir f 2 2 i 2 2\ 

Y = r\i — [r fi + r c — r ) 



2\l/2 



Then 

X = (X + Y) - Y = 2 (rV + r 2 c - r 2 ) 1/2 . (Al) 

We plug the expression for X into Equation [5l noting that, because the shell emits no continuum, 
the maximum value of 9 is not 1, as in the r < r c case, but rather 



/ 2 



2X1/2 



COS 



Thus, 



J(r) 



2( r V + r c 2 -r 2 ) 1/2 dfi, 



where /Uq is the argument of the inverse cosine above. The result is 



J(r) = — < rr c + (r 2 — r 2 ) In 



— \rr c + 



(r 2 - r 2 ) 



In 



^/ r 2 _ r 2 



r + r, 
r — r r 



r + r c 



(A2) 



This preprint was prepared with the AAS IATjjX macros v5.2. 
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The mean intensities inside and outside the core are identical except for two sign differences: 
one in the factor multiplying the logarithm, and the other in the argument of the logarithm itself. 
We find that Eqs. [7] and [8] are continuous at the edge of the core, as they must be, each yielding 
J(r = r c ) = 0.5r c . 



B. Two lines 



We now turn to the two-line case. We denote these lines as R and B, for "red" and "blue." 
First, as stated previously, the source function for the blue line, Sb, is given by Equation [71 The 
source function for the red line, Sr(t), is equal to the mean intensity J\(r), where A is the rest 
wavelength of the line Xo R . The calculation of J is complicated by the fact that photons emitted 
by the blue line may be scattered into resonance with the red line, while some continuum photons 
from the glowing core which would, in the absence of the blue line, redshift into resonance with 
the red line are actually scattered away by the blue line. In velocity space the region where the 
blue line interacts with the red line takes the form of a "scattering sphere" called the common 
point velocity surface (CPVS) (jJeffery fc Branch! I199Q ). The radius of the CPVS, denoted Y in 
Figure \T7\ is given by the Doppler formula, Equation [T4l where one replaces z in that equation 
with Y. In §B.1I we study the interaction of two lines confined to the core, and in §B.2I we explore 
two lines in the shell. 



B.l. Two lines inside the core 



When both lines are confined to the core, Sr contains three components: 



1. the intensity of continuum photons which can redshift into resonance with R, but are scattered 
away when they reach the CPVS; 

2. photons emitted by B along the CPVS which redshift into resonance with R; 

3. continuum photons which form inside the CPVS and therefore redshift into resonance with 
R without interacting with B. 



Segments 1 and 3 are labeled X and Y respectively in Figure [T7J Mathematically we write this as 

Sr = \ j Xe ' TBd ^ + \ J Sb ^ - ^ TB W + \ j Ydfi. (Bl) 

To calculate X = X{(i) we refer to Figure [181 Defining the triangle ABC in Figure [T8l we can use 
the law of cosines to find: 

r 2 c = (X + Yf +r 2 R - 2(X + Y)r R ^ (B2) 
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Fig. 16. — Geometric configuration used to calculate S(r > r c ) in the absence of lines in the core. 




Fig. 17. — The different components of an intensity ray incident on the red line in the two-line 
case. 
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where we have used the fact that X and Y are co-linear and fi = cos9. Solving for X, 



X = (r Rf x -Y)±yJ (run - Y ) 2 + r 2 - r\ + 2Yr Rf x. (B3) 
Expanding the square root yields 



(trim - Y) 2 + r 2 -r 2 + 2Yr R ^ = J r 2 - r 2 (1 - fi 2 ) + Y 2 . (B4) 



When — 1 < /i < Y/r R , we see that (r R fi—Y) < and we must take the positive root in Equation[ 
However, when Y/r R < fi < 1 we see from Equation IB3I that [i > and thus the term in the square 
root is larger than the expression r R fi — Y and so the positive root is also correct. Hence, 



X = (r Rf i-Y) + y/r 2 -r 2 R (l- fJ 2 ) + Y 2 , (B5) 

for — 1 < /j, < 1. Before proceeding we note that if the CPVS extends past the edge of the core, 
that is, if r R + Y > r c , then there exists a critical angle # cr it for which X becomes undefined if 
6 > #crit- Its value is 

Y 2 + r 2 - 1 

Merit = , (B6) 

2Yr R 

where ^ cr ;t = cos C rit- When integrating to find the contribution of the CPVS to the source function 
of the red line, S R , we must stop the integration at this angle. This limit is depicted in Figure [T9l 

Next, we must calculate r$ B = rs B (/J,), the location of the CPVS, as shown in Figure [201 Again 
through vector addition arguments we find 

rsM = (4 + Y 2 ~ ZtrYv) 112 . (B7) 



It is at this location r$ B that both tb and Sb in Equation IB1I are evaluated. Lastly we turn to 
the calculation of Y in Figure [T71 If the entire CPVS fits inside the core, we may use the Doppler 
formula, Equation [111 to calculate Y for all 9. However, if part of the CPVS extends past the edge 
of the core then Y takes a slightly more complicated form. This latter case is shown in Figure [2T| 
where 9\ > 9 CT i t ■ We calculated this angle in Equation IB61 If 9 > 9 CYit then, from Figure EH 

DM = r m ± (r^ 2 + r 2 - r 2 ,) 1 / 2 . (B8) 

Since D is a length and must always satisfy D > we can rule out immediately the "minus" 
solution, since it is negative for all possible values of (r R ,Y,n). Thus we take the positive root: 

D(jj) = r R fi + (r 2 B p 2 + r 2 - r 2 ,) 1 / 2 -!<//< /^it- (B9) 



B.2. Two lines outside the core 



In the case that we have a resonance line that is strong under cold conditions, for example, 
certain lines of Ca II or Mg II, we may want to consider multiple lines forming outside the core. 



Fig. 18. — Explicit geometric construction of two-line configuration used to calculate X = X(fi). 
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Fig. 21 — 



Geometric configuration of truncated intensity ray D = D{fj). 
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Fig. 22. — Case I geometry. The value of A corresponds to the magnitude of the vector xib, E to 
wi, C to v&, D to ut, r c to du, and r s to xh. Only E and C contribute to the intensity of the ray 
since A and D lie outside the emitting core. 
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Figure [22] illustrates the most complex case. If 9 > # C ritj where /^ cr it = cos(# cr it) = 1 — {r c /r R ) 2 , then 
the characteristic in that direction does not intersect the core and accumulates no intensity, but for 
the case that [i > /i cr it then the characteristic is given as shown in Cases I— III, (Figures !22H24p , 

For Case I, Figure [22j we find 



D = r R cos9- Jr^-r 2 R sm 2 9, (BIO) 



C = Y -D = Y -r R cos9 + Jr* - r% sin 2 9, (Bll) 



E = 2(r R cos 9 - D) - C = Jr 2 - r\ sin 2 9 - Y + r R cos 9, (B12) 

and 



A= Jr 2 -r R sin 2 9- Jr 2 -r R sin 2 9. (B13) 
For Case II, Figure [23l we have 



E = 2Jr 2 -r 2 sin 2 6» (B14) 



C = r R cos9 -Y - Jr 2 -r 2 sin 2 9 (B15) 



and 



For Case III, Figure [2 



A = Jr 2 -r R sin 2 9- Jr 2 -r R sin 2 9 (B16) 
D = Y. (B17) 



D = r R cos9- yV 2 - r\ sin 2 9 (B18) 

> (B19) 



'2 — ^.2^ p ir1 ^ ' 



E = Y-D-C = Y-^r 2 - r 2 R sm 2 9 - r R cos9. (B20) 
and 



A = yjr 2 - r 2 R sin 2 9 + r R cos9- Y. (B21) 

Then J for Case I is: 



1 ( /"Merit /-i 

3 = 9 1 / - e " TS ) ^ + / i Ee ~ TB + ^ - e " TS ) + c] ^ r 

^ W~l •'Merit J 



For Case II, J is: 

J f /"Merit /"I 



7 1 r r crit 5b(1 _ e - TB) d/l + r 1 ^ + 5b(1 _ e - TB ^ d l 

1 W-l -'Merit J 



(B22) 



J = ~ I y^ 14 5b(1 - e~ TB ) dp + y [Ee~ TB + S B (1 - e~ rs )] d^ . (B23) 
Finally, for Case III, 



(B24) 



The extension of the two-line case to the iV-line case is straightforward and is implemented 
most easily in a recursive fashion, moving along wavelength as the value of r R increases. 



Fig. 23. — Case II geometry. The value of A corresponds to the magnitude of the vector xw, E to 
wu, C to uv, D to ut, r c to era, and r s to xo. Only E contributes to the intensity of the ray since 
A, C, and D lie outside the emitting core. 



Fig. 24. — Case III geometry. The value of A corresponds to the magnitude of the vector xi, E to 
mi, C to vr&, D to ut, r c to ou, and r s to xo. Only C contributes to the intensity of the ray since 
A, E, and D lie outside the emitting core. 
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B.3. Calculation of emergent intensity 

Given that the flux integral in Equation U] can be written as an integral over impact parameter 
p, we concern ourselves here with the calculation of the emergent intensity of rays with constant p, 
denoted I\(p). We first calculate the domain [z m { n , z max ] over which a ray with given p exists: 

^■min — {j"c P ) ^ 
W = +{rl-p 2 ) 1/2 

We then use the Doppler formula to establish the wavelength domain [A m j n , A max ] over which it is 
possible for a particular wavelength point A to be in resonance in the core, given p: 

Amin = ^ r (B25) 

1 J_ <v I "core \ x 7 

J- T" ^max \ c J 

Amax = VT~y A ftev (B26) 

All lines with rest wavelength Ao which lie in this domain will be in resonance with the ray I\{p) 
at some location in the core. 

Now consider a ray with two lines which both fall in the domain [A m i n , A max ]. The ray redshifts 
into resonance with the blue line "before" (at larger z, closer to the back of the core) redshifting 
into resonance with the red line; we denote these locations zb and zr, respectively, where zb > zr. 
The segment of the ray between the back edge of the core and zb has intensity (r 2 — p 2 ) 1 / 2 — zr- 
However it will redshift into resonance with both lines before emerging from the core, so it acquires 
two attenuation terms, e~ TB and e~ TR , where tb and tr are evaluated at zb and zr, respectively. 
The continuum segment between zr and zr has length zb — zr and is attenuated only by the 
red line. The front-most continuum piece, between zr and the front of the core, is unaffected by 
scattering and has intensity (r 2 — p 2 ) 1 / 2 + zr. 

The source functions of the red and blue lines also contribute to the total emergent intensity. 
The contribution from each is Si(l — e~ n ), where i G {R, B}. However the blue line's contribution 
will be attenuated when it redshifts into resonance with the red line, and so it is receives the usual 
e -Tfl multiplicative factor. Therefore the total emergent intensity is 

hip) = ({r 2 c ~ P 2 ) 1/2 ~ z B )e- TB e~ TR + S B {1 - e - TB )e~ TR + (z B - ZR)e~ TR 

+S R (1 - e~ TR ) + (r 2 c - p 2 ) 1 ' 2 + zr. (B27) 

If for a particular ray one of the two lines is outside [A m i n ,A max ], the result for I\{p) reduces to 
the one- line form for I\(p) of Equation [10] or Eqs. [T5l [T7| and 1181 for lines outside the core. It is 
important to understand that if, for some impact parameter p, the blue line has wavelength outside 
[A m in,A max ], such that I\{p) is given by Equation [TUJ, the source function S = Sr evaluated in 
that equation may nevertheless contain scattering effects of the blue line, as long as at least some 
portion of the CPVS for the blue line lines within the core. 

If both lines fall outside this domain the emergent intensity becomes the pure continuum result. 



